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We propose a dynamic Kerr effect experiment for the distinction between dynamic hetero- 
geneous and homogeneous relaxation in glassy systems. The possibility of this distinction 
is due to the inherent nonlinearity of the Kerr effect signal. We model the slow reorienta- 
tional molecular motion in supercooled liquids in terms of non-inertial rotational diffusion. 
The Kerr effect response, consisting of two terms, is calculated for heterogeneous and for 
homogeneous variants of the stochastic model. It turns out that the experiment is able to 
distinguish between the two scenarios. We furthermore show that exchange between rel- 
atively 'slow' and 'fast' environments does not affect the possibility of frequency-selective 
modifications. It is demonstrated how information about changes in the width of the relax- 
ation time distribution can be obtained from experimental results. 



I. Introduction 

Correlation functions in disordered systems usually decay nonexponentially in time jl]. 
The microscopic origin of this behavior is of importance for a detailed understanding of the 
dynamics in these systems. In principle, the stretched exponential form of the correlation 
functions can be explained by two extreme scenarios. Either different Debye-type relaxation 
rates superimpose to the total response, or the broadened form is an intrinsic property of 
the dynamics in glassy systems. These different scenarios have been denoted as dynamic 
heterogeneous and homogeneous [2]. Several, but not numerous, experimental techniques 
have been developed that allow the investigation of dynamic heterogeneity, like four-time 
NMR measurements [3], a deep bleach experiment ^ and nonresonant hole burning (NHB) 

In a theoretical investigation of NHB [H] it has been shown that indeed heterogeneous 
and homogeneous dynamics can be distinguished. Experimental realizations of NHB have 
been conducted on numerous systems like supercooled liquids |Sl 13 E] , relaxor materials 

and spin glasses [TOj. In all these studies indications for dynamic heterogeneity surpass 
those for homogeneity, but also homogeneous behavior was found, e.g., in an amorphous 
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ion conductor ^T]. Recently a mechanical variant of NHB has been developed |12j . 

Mainly accountable for the distinction between heterogeneous and homogeneous dy- 
namics is the possibility to adress certain dynamics in the sample, i.e., the selection of 
specific dynamic subensembles. Frequency-selective behavior is obtained only if the mea- 
sured response function is nonlinear in the applied field, which holds for the dynamic Kerr 
effect ^2] as well as for NHB. In Kerr effect relaxation, frequency-selective behavior can 
be achieved if a driving AC field with varying frequency is applied. We demonstrated 
frequency-selectivity in a Kerr effect experiment in the terahertz range in a recent pub- 
lication ^3]. Treating the vibrational dynamics around the boson peak in a Brownian 
oscillator model we concluded that the frequency dependence of the vibrations' damping 
can be extracted experimentally. 

The Kerr effect has to our knowledge not yet been used for an investigation of dynamic 
heterogeneities. In this article, we study the Kerr effect response in the range of slow 
reorientational motions in supercooled liquids. We propose the Kerr effect for a distinction 
between dynamic heterogeneous and homogeneous relaxation. The advantage of the Kerr 
effect compared to NHB lies in the fact that the Kerr effect response is purely nonlinear 
in the applied field. Therefore no separation of nonlinear contributions from the linear 
response background is required, and field strengths remarkably weaker than those usually 
applied in NHB (100 kV/cm [71 |H1) niay suffice. It should also be possible to study 
dynamics at different temperatures ranging from the glass transition temperature to high 
temperatures. Here, possible changes in a heterogeneous distribution with temperature 
[T^l Hn] might be investigated. 

The paper is organized as follows. Section II outlines the principles of the Kerr effect, 
explains our suggested experiment and defines heterogeneous and homogeneous models. 
In Section III we give expressions for the Kerr effect response in the different approaches 
and we discuss the results. It turns out that indeed the distinction between dynamic 
homogeneity and heterogeneity is possible. The conclusions are given in Section IV. 

II. Theory 

1. Dynamic Kerr effect 

If an anisotropically polarizable sample is exposed to a time-dependent electric field it 
becomes birefringent. This phenomenon is known as dynamic Kerr effect. In a theoretical 
description the coupling of external electric fields E{t) to matter via the sample's permanent 
dipole moment fi has to be taken into account as well as the coupling via the polarizability 
a. The structure of the Hamiltonian describing the interaction is thus of the form ^Tj 

Hext = -fiE{t)Pi{cose) - aE'^{t)P2icos9) (1) 
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The permanent dipole moment interaction is linear in E, while the polarizability interaction 
is of order E"^ since the induced dipole moment itself is cx: E. The appearance of the 
Legendre polynomials Pl (cos 6) is due to the tensorial nature of the dipole moments. Since 
the permanent dipole moment is a vector and the polarizability is a matrix, they transform 
like first and second rank Legendre polynomials, respectively. The scalar a in Eq.(P3) is to 
be identified with the difference between the polarizabilities parallel and perpendicular to 
the internal symmetry axis [18], and 6 is the angle between this axis and the applied field. 

The time-dependent polarizability is the quantity of interest for a description of the 
dynamic Kerr effect. We therefore have to calculate the expectation value of the second 
Legendre polynomial 

{aP2{cose{t))) =a{P2{t)). (2) 

The brackets denote an expectation value over the whole sample. In this work we focus on 
reorientational motions. Here, the time dependence of the orientation (described by the 
angle 6) determines the time dependence of the polarizability. In the following, we will use 
the short hand notation on the right hand side of Eq.Q. 

The expectation value {P2(t)) is calculated in some dynamic model, where the applied 
external field that determines the time dependence is treated in perturbation theory. Even 
without specifying the model yet, it is clear that the linear response of order E must vanish 
as long as isotropic systems are considered. This is because the linear response is always 
proportional to the first rank Legendre polynomial, see Eq.(^. Because of the orthogonality 
of the Legendre polynomials the linear response vanishes if the expectation value of the 
second rank Legendre polynomial, Eq.(j21), is calculated. The Kerr effect response is thus 
of order E"^. 

In some more detail we have two contributions to the expectation value {P2{t)). The 
first one is proportional to the polarizability a (formally calculated with the second term in 
the external Hamiltonian in first order perturbation theory), the second one is quadratic 
in the permanent dipole moment fi (second order perturbation theory with the first term 
in (^). We will denote these contributions in the following as {P2{t)) and {P2^{t)). 



2. Experiment 

We propose an experiment as follows. First, a sinusoidal electric field is applied for an 
arbitrary number of (full) cycles to a sample in equilibrium. 

2'kN 

E{t) = E sm{%T) T<tp = —— (3) 

i Lp 

The most important experimental parameter is the pump frequency Qp. After the 'pump 
time' tp the field is switched off and the relaxation of the sample's polarizability back to 
thermal equilibrium is measured. In our terminology the time variable t starts with at 
tp. 
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The experiment is based on the same idea as NHB. If the dynamics is assumed to 
be heterogeneous (different relaxation rates), then it should be possible to address 'slow' 
or 'fast' dynamics separately by varying the pump frequency from 'low' to 'high'. The 
so-selected ensembles will then show 'slow' or 'fast' relaxation behavior afterwards. In 
a Fourier transformed representation of the relaxation, one would therefore expect an 
extremum position proportional to the applied field frequency. On the other hand, in a 
homogeneous system one would expect relaxation on a single time scale only, and different 
pump frequencies should not alter the time-dependence or extremum position in the Fourier 
transform. Thus, shifts in the extremum position can be considered as indications for 
dynamic heterogeneous relaxation while no shifts would support the homogeneous picture. 

The aim of this paper is to demonstrate that the concept outlined above indeed holds in 
heterogeneous and homogeneous models. Because of the nonlinearity of the Kerr effect the 
contributions {P2{t)) and {P2^{t)) are not trivially related to the correlation function, and 
have thus to be explicitly calculated for a given field sequence like Eq.Q. We specify the 
dynamic models used in the calculation of the response functions in the following section. 



3. Dynamics in the rotational diffusion model 



For a description of the dynamics we use the well-known rotational diffusion model [TH] . 
The rotational diffusion equation for the conditional probability Pi{Q,t\Q' ,t') to find the 
orientation Q at time t, assuming it was Q' at time t' < t is 



d_ 

dt 



p^{n, t\n', t') = iii{n, t)Pi{n, t\n', t') 



with the Fokker-Planck operator 
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Here, 9 and denote polar and azimuth angles with respect to the molecules internal 
symmetry axis. These angles constitute the orientation VL. Di{t) is a (possibly time- 
dependent) rotational diffusion constant, (3 = {ksT)'^ and Tiext is the Hamiltonian for the 
interaction with the applied external field in Eq.(|T|. 

The solution of Eq.© without driving field {Tiext = 0) is achieved via an expansion in 
Wigner rotation matrices [20] • The ansatz 



Pi{n,t\n',t') = ^ 

inserted in Eq.(jll) yields 



LMN 



(6) 



GP(t,t')=exp 



-L(L + 1) f'dr D,{t) 
Jt' 



(7) 
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and rotational correlation functions of the L-tli Legendre polynomial are given by 

C^{t,t') = {PL{cose{t))PL{cos9{t'))), = ^±-^G\'^\t,t'). (8) 

In the presence of external fields perturbation theory is used to calculate the response func- 
tions. The corresponding perturbation operators result from the Fokker-Planck operator 
(0) including the external Hamiltonian in Eq.([Q). Linear response functions obtained this 
way are related to the correlation function via the fluctuation dissipation theorem pij . 



Heterogeneous model 

In a heterogeneous scenario one assumes the coexistence of different, 'slow' and 'fast' 
environments related to small and large rotational diffusion constants Di. If Di is time- 
independent, the propagator Pi{Q,t\Q' ,t') is time-translational invariant and the correla- 
tion function in case of a single rotational diffusion constant is given by 

C^{t,t') = C^it - t') = ^^exp [-L{L + l)A(t - t')] . (9) 

Dynamic heterogeneities are known to have finite life times of similar magnitude as the 
a-relaxation time fl^ EB] ■ Finite life times can be modelled via exchange models. Various 
approaches exist to include exchange in a model |21] . In a master equation ansatz |2S1 one 
assumes that slow environments become fast and vice versa. The environments are related 
to rotational diffusion constants, and exchange between the environments leads to effec- 
tively time- dependent rotational diffusion constants. Time-translational invariance is pre- 
served within this framework. We will denote the environments by in the following, and 
we add a second index to quantities like the propagator. Pij{Q,t\Q' ,t') = Pij{Q,t — t'\Q' ,0) 
thus is the probability to find the orientation Q in environment at time t, assuming it was 
Q' in environment ej at time t' . Transition rates jij determine transitions from environment 
ej to ej, meaning that the rotational diffusion constant changes its value from Dj = D{ej) 
to Di = D{ei). Furthermore, the behavior of the orientation during a transition has to be 
specified. If we assume that each transition is associated with a rotational random jump, 
then an exchange model can be defined by [23] 

^Pki{n,t\Q',t') = (n,(fi) -7,) P,,{Q,t\Q',t') + -L^7,^. f dQoPA^o,m',t'). (10) 
at m .^^ J 

Here, nfc(n) is the Fokker-Planck operator (0) with Di{t) replaced by the time- independent 
and = J2jjtkljk is the total rate out of state e^. 

The solution of Eq. pO|) without external fields is obtained in a similar manner as the 
solution of Eq.(jD). Here, the function in Eq.® is replaced Gi^\t,t'). For L 

the elements of the matrix G^^\t,t') are given by 

Glf°)(t,t') = feexp[AP(t-t')l (11) 



5 



with the eigenvalues 

Af) = -L(L + l)A-7i- (12) 

The equihbrium probabihties Pf^ to find the system in environment in the pres- 
ence of exchange (finite jij) are obtained from the long-time limit of the propagator, 
P/'' = \im(^t-t')^oo I d^Pik{^,t\^l' ,t') = gI^IoOjO), which guarantees detailed balance. 
Correlation functions in the heterogeneous model are then calculated as 

c'{t) = Y.cHt)pr = ^^E^xp [A^h] pr. (13) 

i ' i 

The larger the exchange rates are, the more they dominate the eigenvalues in comparison 
to the L-dependent term in Eq. fjl2p . This also holds in other exchange models [21]. Large 
exchange rates thus lead to a weaker L-dependence of the eigenvalues, and the exponents 
in correlation functions ()13p for different values of L become more and more similar. The 
inclusion of exchange in the rotational diffusion model thus leads to very similar correlation 
times for different L, as usually found experimentally [22] • The calculation of response 
functions is straightforward. 

Exchange effects are neglected if 'jij = is chosen. Then the eigenvalues in Eq. ()13|l 
are given by A^-^^ = —L{L + l)Di, and a distribution of rotational diffusion constants for 
a heterogeneous model can be arbitrarily chosen. If the relaxation times Xj oc D^^ are 
distributed according to the generalized Gamma distribution p!U], a stretched exponential 
form of the correlation function is obtained. 



Homogeneous model 

In a homogeneous scenario one assumes that the relaxation rate (or rotational diffusion 
constant) depends on time. A decaying power law function 



It 

inserted into the correlation function (jH)) using Eq.© yields 



(14) 



(Pi(t)Pi(t')) = -exp {-{^Irf- 



exp 



(15) 



Thus, the special choice in Eq. (ll4j) reproduces a stretched exponential for the correlation 
function. Note that Eq.(|15p is not time-translational invariant, which is due to the time 
dependence of -Di(t) j^. Due to the lack of time-translational invariance the system has 
an 'internal clock' that is hard to justify from a physical point of view. The consequences in 
our case will be discussed below. An alternative for modelling homogeneous dynamics are 
the so-called fractional Fokker-Planck approaches ^T]. Here, time-translational invariance 
is maintained 
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III. Results 



In this section we discuss the Kerr effect response to the alternating field in the various 
approaches considered in the previous section. After the orientational integrations have 
been carried out we find the expressions 

(p-(t)) = -[3aE^ y dt' sm\nj)D,{t')Gf{t + t„, t')PP (16) 
5 

for the aa- contribution to the Kerr effect response and 

5 



X ^ 
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V dt' t dt" sin(n/) sin{npt")Di{t')D,{t")Gf\t + t„, t')Gf\t\ t")Pr (17) 
, Jo Jo 



for the a/i/i-contribution. Eqns. p6ll7p are valid for heterogeneous and homogeneous mod- 
els. In the heterogeneous case the Di are time-independent and in the homogeneous case 
no summation over i is required. In the presence of exchange G^j^\ti^t2) is to be identified 
with the diagonal elements of the matrix G*-''"^(ti, ^2), see Eq.()ll|). 

Note that Eq.(|T^ is formally identical to the expression for the linear (dielectric) re- 
sponse, except for the fact that the external field of sine form enters quadratically and 
Gf^ appears instead of Gf^ . As a consequence of this quasi-linearity the aa-contribution 
is -like the dielectric linear response- unable to distinguish between homogeneous and het- 
erogeneous scenarios. 

Let us briefly discuss the simplest case where the system is described by a single, time- 
independent rotational diffusion constant. Here, G^j^\ti,t2) are exponential functions ac- 
cording to Eq.Q. We furthermore have to set Di = const and J2i Pi'^ = 1 in Eqns. ()lfipi7j) . 
The time integrals can then be carried out analytically, leading to response functions of 
the form 

{Pf{t)) = v4®(fip. A, tp) exp(-6At) (18) 

where B stands for a or for yU/i. In the following discussion of the response functions we will 
show the imaginary part of Fourier-transformed response functions, denoted as {Pf{uj)). 
The symbol Umax will be used for the extremum positions in {Pf{uj)). 

Both contributions ()18j) with Q = a and G = /i/i show the same time dependence, a 
monoexponential decay with the decay rate of the second Legendre polynomial. This is 
analogous to the step field-off- response discussed in Ref.[2^. The maximum position in 
{Pf{(jj)) therefore is located at frequency uJmax = S-Dj for both terms. 

The amplitude functions A^^flp, Di,tp) are given by jHOj 

A°(fip,A,tp) = ^PaE^^e"{6D,,2np)[l-e~'''''-] (19) 
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X (3e"{4Di, Qp) [e'^-^^*'' - e"'^'*''] + -e"{6Di, 2Qp) [l - g-*^^'*^]) (20) 

with the dielectric loss defined as e"{x,y) = xyjix^ + y^). In the context of frequency- 
selectivity the behavior of the amplitudes as a function of pump frequency is important. 
A^^i^p^ Di, tp) plotted versus pump frequency shows a peak because the amplitudes become 
large if the arguments in the dielectric loss functions e"{x. y) are equal, that means for 
VLp 2Di, 3Di, 4:Di. If for a moment we neglect the fact that the arguments of the dielectric 
loss functions have varying prefactors, then E'^e"{Di, flp^ represents the pump frequency- 
dependent energy absorbed by the system. 

No frequency-selective behavior can be expected from the quasilinear aa-term, where 
the additional factor Jlp appearing in the nominator of A'^{D,p, Di,tp) does not allow for 
a 'resonant' energy absorption. This is because the product Qpe"{6Di, 2Qp) is a strictly 
increasing function of Jlp in contrast to e"{6Di, 2D,p) which exhibits a peak at ftp — SDi. 

Heterogeneous model 

In a first discussion of Kerr effect responses in the heterogeneous model we neglect exchange 
effects. We choose the rotational diffusion constants distributed according to a general- 
ized Gamma distribution with parameters such that the resulting (normalized) (L = 1)- 

-(t/r)^^j with the KWW-exponent f3K = 0.5 and r = 1. 
diffusion constants is shown in the upper panel of Figure 
1. The lower panel shows the Kerr efiiect signal {P2^(uj)) for different pump frequencies 
Q,p — 0.01, 0.1, 1, 10, 100 (indicated by letters A,B,C,D,E in the figure). One cycle has been 
applied to the system for all pump frequencies. 

In order to understand the curves in Figure 1 we refer to the previous section where 
we discussed the response of a single Di. Here, the energy absorption was shown to have a 
maximum for Vlp ^ 2Di,3Di,ADi. This means that in a distribution of Di most energy is 
absorbed if Di ^ Qp/2 holds. Thus, a pump frequency dependent selection is possible. The 
Di corresponding to Qp = 10^^ . . . 10^ are marked by the vertical arrows in the distribution 
function (upper panel of Fig.l). Under the assumption that the response is mostly governed 
by the selected Di we obtain a maximum position at ujmax — ^Di, thus at comax ~ 3Qp. 
Although the prefactor 3 is not exactly recovered in the extremum positions in the lower 
panel of Fig.l, we indeed observe a proportionality of extremum position (indicated by 
dotted lines) and pump frequency for flp — 100... 0.1 (E,D,C,B). These shifts of the 
extremum position are the hallmark of dynamic heterogeneity. For Qp = 0.01 (A) the 
relative shift becomes weaker. Here, a selection of Di ^ 0.01/2 does not determine the 
total response because these small Di are underrepresented in the distribution function, 
see the arrow (A) in the upper panel. In that case, rotational diffusion constants around 



correlation has the form exp 
The distribution of rotationa' 
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Di pa 0.05 (B in the upper panel) determine the extremum position. This also holds for 
still smaller pump frequencies. Here, the extremum position becomes independent of VLp 
while the amplitudes of the signal strongly decrease. A similar behavior is observed at the 
high-frequency end of the distribution, which we do not show explicitly. The maximum 
position as a function of Vtp will be discussed in detail later compared to results in the 
homogeneous limit. 

In order to discuss the width of the responses, we added a Debye function (dotted line 
near C in the lower panel of Figure 1). The heterogeneous response (C) is somewhat broad- 
ened in comparison. The curve for Vlp = 100 (E) has a negative onset on the low-frequency 
side of the peak. Two remarks concerning this finding of a partially negative response 
function are appropriate. First, the negative sign of the response function is due to the 
first term in Eq. ()20|) . which is proportional to the difference [exp(— 6Djtp) — exp(— 2Djtp)]. 
This is negative for finite tp in contrast to the other term oc [1 — exp{—6Ditp)]. A negative 
response function therefore is a switch-on effect, since the former contribution vanishes for 
large tp. Indeed, if several cycles are applied, all response functions become positive in the 
whole frequency range. Second, even for finite tp, the negative term can be considered as 
an artefact of the rotational diffusion model. This is because the rates QDi, 2Di in the 
exponents are the decay rates of the second and the first Legendre polynomial, which are 
known to be of similar magnitude in real systems as discussed in Section II. In any model 
where these rates are very similar (e.g. the exchange model discussed in Section II) or 
equal (e.g. isotropic rotational jumps instead of rotational diffusion) no negative sign of 
the response occurs for arbitrary tp. 

Figure 2 shows the contributions of the quasilinear term (P|*(c(j)) with the same pa- 
rameters as in Figure 1. The extremum position is not proportional to pump frequency 
although it weakly depends on Qp. While the pump frequency spreads over five decades, 
the maximum position shifts only for less than one decade. The shift of the extremum po- 
sitions is further diminished if several cycles are applied. This demonstrates the expected 
result that the quasilinear term is not able to detect djTiamic heterogeneities. In the fol- 
lowing discussion, we will therefore focus on the a/i/i-contribution. The quasilinear term 
can be viewed as a background contribution to the response function. As it is of quadratic 
order in the applied field like the a/z/z-contribution, no phase-cycling procedure allows the 
extraction of either of the contributions. However, the terms should be easily distinguish- 
able since the quasilinear signal is at a more or less constant position. Furthermore, the 
problem with the quasilinear term is trivially minimized if systems with large permanent 
dipole moment /i (compared to their polarizabihty a) are studied. In this case, (P|'(^^)) is 
much smaller than {P2^{uj)). 

We now discuss the influence of exchange on the Kerr effect response in the frame- 
work of the exchange model outlined in Section II. For our calculations we assume that 
the exchange rates factor into terms that depend on initial and final environment only, 
7jj = p(ej)X exp(/3ej). Here, p(ej) (probability to end in final state e^) is some normalized 



9 



density of states and Xexp(/3ej) (probability to leave initial state e^) can be interpreted as 
thermally activated jumps out of environment ej. Here, we have identified the so-far not 
specified ej with energies. The prefactor X allows to adjust the 'strength' of exchange. We 
furthermore relate the rotational diffusion constants to as Di = Yexp{i3ei). The model 
is designed such that the total rate out of environment is proportional to the appropriate 
rotational diffusion constant, 7j ^ {X/Y)Di, meaning that slow environments exchange 
slowly and fast environments exchange fast. The ratio of 7^ and Di, given by {X/Y), 
determines the influence of exchange on the decay of rotational correlation functions. 

For the curves shown in Figure 3 we chose p(ej) to be Gaussian with width a = 3 and 
we have set (3 = 1. In order to obtain correlation times around unity we set Y = 100. 
The upper panel of Figure 3 shows normalized correlation functions for L = 1,2 without 
exchange [X/Y = 0, dashed lines) and in the presence of exchange {X/Y = 10, solid lines). 
In the latter case the decay is strongly affected by exchange. All correlation functions shown 
can be well fitted with a function of KWW-type with a stretching exponent ~ 0.5 . 
The ratios of correlation times obtained from the fit are ti/t2 ~ 3.0 without exchange 
and T1/T2 ~ 1.3 with exchange. 

The nonlinear Kerr effect response {P2^{ijj)) in the presence of exchange is shown in 
the lower panel of Figure 3. All parameters describing exchange are chosen as for the solid 
lines in the upper panel. In the limit without exchange, {X/Y) = 0, the model discussed 
in Figure 1 is analytically recovered with a modified distribution of rotational diffusion 
constants. The resulting curves in that case are very similar to the responses in Fig.l. For 
that reason we show only the curves with exchange {X/Y = 10). The pump frequencies 
are chosen as Qp = 10^^ . . . 10^ with one applied cycle for each frequency. The shifts 
of the extremum position obtained are very similar as in the heterogeneous case without 
exchange, cf. Fig.l. The peaks become slightly broader in the presence of exchange. No 
negative sign of the response occurs in any of the curves. 

We therefore conclude that exchange does not affect the possibility to detect dynamic 
heterogeneities. It is important to stress that the Kerr effect experiment does not allow to 
conduct information regarding the lifetime of dynamic heterogeneities. 



Homogeneous model 



We now discuss the responses in a homogeneous scenario, where the distribution of rota- 
tional diffusion constants is replaced by a single, but time-dependent rotational diffusion 
constant. The time dependence is chosen in a manner that the correlation function for 
L = 1 is proportional to exp —{t/rY^ with (3k = 0.5 and r = 1. The correlation func- 
tion is thus the same as in the situation we discussed in Figure 1. Inserting Eqns.(|7jl. (|T^ 
into (fTTjl we find 



exp 



-3 



t + t. 



T 



X 



10 



1 








exp 




exp 













.(21) 



The double integral does not depend on time t and therefore represents only a prefactor to 



the response. The time-dependence is solely governed by the function exp 

Thus, the maximum position in {P2^{ijj)) is determined by the Fourier-transform of this 
function. 

Figure 4 shows the homogeneous responses for the same pump frequencies as in Figure 
1. One cycle has been applied for each pump frequency, tp = 27r/Qp. It is this pump 
frequency-dependence of tp which leads to the shifts in the extremum positions in Figure 4 
(see the function above). Although the cxtrcma arc not at a constant position as usually 
expected in a homogeneous scenario, the shifts are less pronounced than in the heterogenous 
limit. The dashed curves corresponding to flp — 1 and flp = 10 are plotted negative. This 
is an example for the fact that signs of nonlinear responses are in general not easy to 
predict. The change of sign is not a systematic feature in our calculations. 



Comparison of heterogeneous and homogeneous models 

In Figure 5 we compare the maximum position in heterogeneous (upper straight line) 
and homogeneous (lower straight line) models, plotted versus pump frequency Qp for one 
applied cycle. The distribution in the heterogeneous case is as in Fig.l, and the time- 
dependent rotational diffusion constant of the homogeneous model is as in Fig.4. In the 
homogeneous case a maximum position independent of pump frequency is obtained only 
for sufficiently large Qp. Clear shifts of the extremum position arc observed in the hetero- 
geneous limit. The shifting becomes weaker for Qp < 0.1. This is due to the shape of the 
chosen distribution function as elucidated in the discussion of Fig.l. 

For a distinction between homogeneity and heterogeneity, the straight lines in Figure 
5 do not look too promising at first glance because the homogeneous extrema are not 
at a constant position. However, this finding should not be overemphasized. The ob- 
servation of shifts in the homogeneous extremum positions relies directly on the lack of 
time-translational invariance of the specific homogeneous model used in the calculation. 



-3 



t+tp 
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The fact that the 



This leads to the shifted time (t -|- tp) in the function exp 

extrema arc not at a constant position is only due to the appearance of the pump time tp 
in this function, which is related to the pump frequency via tp = 211 /Q.p. It is furthermore 
easy to eliminate this effect. If tp is chosen independent of the pump frequency, then the 
homogeneous model obviously yields extremum positions completely independent of the 
pump frequency. This can be achieved by applying several cycles N{D,p) ex because of 
tp = 2T:N{Vtp)/Vtp. 

We added some extremum positions for that case in Fig.5, denoted by symbols. We 
chose tp — 27r/0.01, meaning that one cycle is applied for the smallest pump frequency 
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considered, flp = 0.01, and 10^ cycles are applied for the largest Qp = 100. Then the 
homogeneous model extrema (open symbols) do not depend on Qp. The influence of several 
applied cycles on the extremum positions in the heterogeneous case is small. The extrema 
are located at somewhat smaller frequencies in comparison to the curve for one applied 
cycle. This is due to the pump time dependence of the function A^^f^{Qp, Di,tp). 

The theoretical finding that the homogeneous ansatz leads to strictly constant ex- 
tremum positions only if the pump time is constant for all applied frequencies also holds 
for NHB. If NHB is carried out as outlined in Ref. [HT] and the dynamics are treated in 
a rotational diffusion ansatz as in this work, we find in the homogeneous limit the same 
expression as Eq. (j2ip up to prefactors and a modified time dependence. It is again the 
lack of time-translational invariance that leads to terms like (t + tp)^^ in the appearing 
exponential functions. We therefore suggest also in NHB to apply several cycles for higher 
pump frequencies. Also the heterogeneous model for NHB leads to very similar results 
as presented in the previous sections. Shifts of the extremum position are observed if the 
pump frequencies are varied in the range defined by the distribution of rotational diffusion 
constants. 

Distribution width in the heterogeneous model 

Finally, we discuss the influence of the width of a distribution of Di in a heterogeneous 
model. A linear behavior of the extremum position as a function of pump frequencies, 
(^max fl^ with 7] = 1, is obtained only if the distribution of rotational diffusion constants 
is sufficiently fiat in the studied range. If the distribution is narrow, then usually exponents 
1] < 1 are found. 

In Figure 6, extremum positions obtained from the heterogeneous model with 
a logarithmic Gaussian distribution of rotational diffusion constants g{Di) oc 
exp (— [InDj — lnDoo]^/[2cr^]) are shown for three different distribution widths a = 1.8 
(diamonds), a = 1.2 (squares) and a = 0.6 (circles). Exchange effects are neglected. In 
each case, Doo was chosen such that the correlation time of the {L = l)-correlation function 
in a KWW-fit is ri = 1. The corresponding KWW-exponents obtained from the fits of 
the correlation function are Pk = 0.51 (diamonds), Pk = 0.67 (squares) and Pk = 0.88 
(circles). We applied the field for fixed pump time tp = 27r/0.001 independent of the pump 
frequency. 

Heterogeneities are observed via shifts of the extremum position. Constant behavior 
of uJmaxi^p) is obscrvcd for small and large Qp. This is because, similar as in Figure 5, 
the 'selected' in these cases are located beyond the main contributions to the distribu- 
tion functions. In Fig. 5 a constant behavior for large flp is not reached in the frequency 
range considered due to the tailing of the corresponding distribution of rotational diffusion 
constants, see the upper panel of Fig.l. We do not expect that a constant extremum po- 
sition for large Qp would be obtained experimentally. Here, other dynamics than isotropic 
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reorientations will occur in real systems. 

One can see that the rise of 0J„iax with Vtp is the steeper the broader the distribution 
is, approaching the limit 77 = 1 for a very flat distribution. Acceptable power law fits for 
the shown curves are obtained in the window Q,p — Q.l... 10. Here, we find the exponents 
77 f« 0.84 (diamonds), 77 0.63 (squares) and 77 ^ 0.2 (circles). 

IV. Conclusions 

We have calculated the Kerr after-effect response to an external oscillating electric field 
for reorient at ional dynamics in an isotropic system. Such an experiment is very similar to 
NHB, thus aiming at the question of the nature of the dynamics. The different assumptions 
of heterogeneous and homogeneous dynamics have been treated in a rotational diffusion 
model. We find that homogeneous and heterogeneous dynamics show different behaviour 
in the Kerr effect response function. 

We propose the Kerr effect as an alternative to NHB. NHB and Kerr effect have in 
common that they are both nonlinear experiments. The advantage of the Kerr effect is 
that the signal has no linear response background. In an experiment, one would measure 
the reequilibration of the polarizability after the sample has been excited by a sinusoidal 
electric field. The frequency of the latter is then varied. If the time-dependent functions 
obtained in this way are plotted in a Fourier-transformed representation (imaginary part), 
then extremum positions as a function of the pump frequency shift under the assumption 
of heterogeneous dynamics, and they are constant in the homogeneous limit. The Kerr 
effect signal is thus able to distinguish between dynamic heterogeneity and homogeneity. 

The interesting relation between dynamic and spatial heterogeneities is not resolved 
yet. We therefore emphasize that the Kerr effect experiment outlined in this article does 
not allow to obtain any information about spatial aspects of the heterogeneities. 

The influence of exchange, and thus a flnite lifetime of heterogeneities on the responses 
is small. Exchange in our model does not affect the possibility of frequency selective 
modification. 

It turned out that the general trends of homogeneous and heterogeneous limits are 
obtained more clearly if several cycles of the pump field are applied. This is because switch- 
on effects are minimized for longer pump times. In addition, the extremum positions in our 
homogeneous model depend on the pump time tp = 27iN/Qp, and therefore on the pump 
frequency. In order to suppress this artificial dependence, we suggest to keep the pump 
time constant rather than the number of cycles applied if measurements with different Qp 
are compared. If, e.g., N cycles (where ideally N > 1) are applied for pump frequency flp, 
then 2N cycles should be applied if the pump frequency is 2ftp. 

If the width of a relaxation time distribution decreases then the shifts of the extremum 
position with pump frequency become weaker. This point might be of interest in the 
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investigation of a sample at different temperatures. A steeper rise of the extremum position 
with pump frequency for lower sample temperatures would indicate an increasing width 
of the corresponding distribution function. A point we have not discussed explicitly is 
that of course also combinations of heterogeneous and homogeneous behavior may occur 
in real systems. If the weight of homogeneous and heterogenous character of the dynamics 
changes with temperature this should be observable in more or less pronounced shifts of 
the extrema. 
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Figure captions 

Fig.l : Upper panel: The generalized Gamma distribution plotted as a function of rota- 
tional diffusion constants Di — (2td)~^, where td are the more common relaxation 
times. Parameters are chosen such that the correlation function C^{t) has the form 
exp(-(t/r)^^) with r = 1 and (3k = 0.5 . 

Lower panel: Imaginary part of the Fourier transformed Kerr effect response (contri- 
bution of order a/x/x) in the heterogeneous model with the distribution of rotational 
diffusion constants shown in the upper panel. The pump frequencies are Jlp = 10" 
with n — — 2,— 1,0,1,2 (left to right). One cycle of each frequency has been applied 
to the system. The prefactors P'^fi^E'^ were set to unity. For details see the text. 

Fig. 2 : As in the lower panel of Figure 1 for the Kerr effect contribution of order aa. 
Pump frequencies are — 0.01 . . . 100 (extrema from left to right) with one applied 
cycle. Here, PaE'^ was set to unity. 

Fig.3 : Upper panel: Normalized correlation functions with exchange (solid lines) and 
without exchange (dashed lines) for the values of L = 1 and L = 2 in both cases, 
plotted versus time. The respectively faster decaying curve is for L = 2, the slower 
decaying curve for L = 1. 

Lower panel: Imaginary part of the Fourier transformed Kerr effect response (contri- 
bution of order a/Xyu) in the exchange model, with parameters as for the solid lines in 
the upper panel. Pump frequencies are Qp = 0.01 . . . 100 (extrema from left to right) 
with one applied cycle. Here, P^/i'^E'^ = 1. For details see text. 

Fig.4 : Imaginary part of the Fourier transformed Kerr effect response (contribution 
of order an/j.) in the homogeneous model. Pump frequencies are flp = 0.01 . . . 100 
(extrema from left to right) with one applied cycle. Dotted lines are plotted negative. 
We set (3^ii^E^ = 1. 

Fig. 5 : Extremum positions in the imaginary part of the Fourier transformed Kerr effect 
response of order afi/i, plotted versus pump frequency. The lines are for one applied 

cycle of each pump frequency in heterogeneous (upper solid line) and homogeneous 
cases (lower solid line). Connected symbols show the extremum position in the 
case of fixed pump time tp = 27r/0.01 (full symbols: heterogeneous, open symbols: 
homogeneous). For details see text. 

Fig.6 : Extremum positions in the imaginary part of the Fourier transformed Kerr effect 
response of order afifi, plotted versus pump frequency. The pump time is tp = 
27r/0.001 for all pump frequencies. The three curves represent different widths of 
the distribution of rotational diffusion constants with the values a = 0.6 (circles), 
a — 1.2 (squares) and a — 1.8 (diamonds). For details see text. 
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Haberla^Diezemann, Figure 1 




Haberle, Diezemann, Figure 2 
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Haberl^° Diezemann, Figure 3 




Haberle, Diezemann, Figure 4 
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